Spin current generation from Coulomb- Rashba interaction in semiconductor bilayers 



M. M. Glazov/ M. A. Semina,i S. M. Badalyan,2'3,4 q Vignale^ 

^ loffe Physical- Technical Institute of the Russian Academy of Sciences, St. Petersburg, 194021 Russia 
^School of Physics, Astronomy, and Computational Sciences, 
George Mason University, Fairfax Virginia 22030, USA 
^Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 
^Department of Physics, University of Regensburg, 93040 Regensburg, Germany 
^Department of Physics and Astronomy, University of Missouri, Columbia, Missouri 65211, USA 

Electrons in double-layer semiconductor heterostructures experience a special type of spin-orbit 
interaction which arises in each layer from the perpendicular component of the Coulomb electric 
field created by electron density fluctuations in the other layer. We show that this interaction, acting 
in combination with the usual spin-orbit interaction, can generate a spin current in one layer when 
a charge current is driven in the other. This effect is symmetry- wise distinct from the spin Hall 
drag. The spin current is not, in general, perpendicular to the drive current. 



Spin-orbit interactions in semiconductors are tradi- 
tionally studied within a one-electron picture, but there 
are instances in which electron-electron interactions latch 
onto spin-orbit interactions to produce intriguing effects, 
which may lead to the creation of radically new spin- 
based electronic devices [1-3 . Of particular interest is 
the generation of spin currents and spin accumulation 
by an electric current, through the so-called spin Hall 
effect [4]. Recently, a special type of spin Hall effect 
was predicted to occur in double-layer heterostructures, 
i.e., two parallel quantum wells separated by an essen- 
tially impenetrable potential barrier, with a quasi- two 
dimensional electron gas in each layer. This effect, called 
"spin Hall drag" [6], consists of the generation of trans- 
verse spin accumulation in one layer by an electric current 
flowing along the other layer and is caused by the com- 
ponent of the Coulomb electric field parallel to the two 
layers. In this Letter we study a different and novel effect, 
which is driven by the component of the Coulomb electric 
field perpendicular to the layers. This perpendicular field 
creates an inhomogeneous Rashba spin-orbit interaction, 
with spatial variation in the plane of the layer [J . We 
refer to this interaction as the "Coulomb-Rashba inter- 
action" . 

The system under study is shown in Fig. [T] A steady 
electric current is driven in the active layer 2. We show 
that the interplay of the Coulomb-Rashba interaction 
with the ordinary cubic Dresselhaus spin-orbit interac- 
tion (characteristic of semiconductors of the zincblende 
structure), provides a new mechanism for the generation 
of a spin current in the passive layer 1. This is partic- 
ularly remarkable in view of the fact that ordinary spin 
Hall eflFect and spin Hall drag are suppressed by spin pre- 
cession. But, in this case the presence of spin precession 
is absolutely essential to the effect. 

The spin current generation can be understood as a 
two-stage process, which is schematically described in 
Fig. [2] [10^. In the first stage, the steady current of elec- 
trons in layer 2 induces a quadrupolar spin distribution 
in layer 1. To understand qualitatively how this comes 
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Figure 1. Scheme of the device under study: an electric field, 
E, is applied to the active layer (2), a spin current is generated 
in the passive layer (1). Blue circles depict electrons, and the 
wavy line shows the inter-layer Coulomb interaction. Arrows 
marked fci, k2 and pi, p2 are the wave vectors in the initial 
and final states. The vertical arrow emphasizes the relevant 
component of the interlayer Coulomb field. 

about, observe that Coulomb collisions between electrons 
in the two layers take place on average with a positive 
momentum transfer from layer 2 to layer 1 along the 
X axis. Moreover, due to the spin-dependent terms in 
the scattering rate (see Eq. [l] below) , the scattering ef- 
ficiency depends on the relative orientations of electron 
spin and its initial and final wavevectors, denoted by k 
and p respectively. Consider electrons with spins parallel 
or antiparallel to the y axis: Sy > and Sy < respec- 
tively. For Sy > the scattering rate is maximal for 
kx -^Px > and minimal for kx -\-Px < 0. For Sy < the 
situation is reversed: the scattering rate is maximal for 
kx -\- Px < ^ and minimal for kx -\- Px > The strongest 
transitions are marked by solid arrows in Fig. [2|a). As 
a result, states with large values of \ky\ {\ky\ ~ kp) be- 
come depleted of 5^ > electrons and filled with Sy < 
electrons. At the same time, states with large values of 
\kx\ {\kx\ ^ kp) are depleted oi Sy < ^ electrons and 
filled with Sy > electrons. Hence, a quadrupolar spin 
distribution is formed, as shown in Fig. [2|b). 

The second stage of the spin current generation is re- 
lated to the spin precession in the cubic Dresselhaus field. 
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Figure 2. Schematic illustration of the spin current generation 
in the passive layer J^. The circle is the Fermi surface, the 
arrows are the electron spins. Panel (a) shows the scattering 
stage of the process: top part shows electrons with spin par- 
allel to the y axis, bottom part shows electrons with spin an- 
tiparallel to the y axis. The arrows show the spin-dependent 
scattering process: solid arrows indicate the stronger transi- 
tions, having kx + Px > for > and kx + Px < ^ for 
< 0, while dashed arrows indicate the weaker transitions. 
Scattering processes that increase the x component of the 
wave vector (i.e. with Qx > 0) dominate due to the current 
flowing in the active layer, (b) Precession stage of the process. 
The resulting spin distribution after the scattering contains 
second angular harmonics. Green arrows demonstrate the 
spin precession caused by the Dresselhaus field fife, Eq. Q. 
Electron spins with opposite wave vectors precess in opposite 
directions. The resulting dipolar distribution of electron spins 
is presented in the panel (c). 



This is also illustrated in Fig. |2|b), where the green ar- 
rows show the direction of the spin precession induced 
by the cubic Dresselhaus field, fi^, (see Eq. |4] below). 
The field tilts the spins out of the x — y plane, thus cre- 
ating a dipolar distribution of the z-component of the 
spin as shown in Fig. |2|c). As a result, a z-spin current 
is formed. The characteristic C2v symmetry of systems 
with both Rashba and Dresselhaus interactions causes 
the spin current to be parallel to the driving electric field 
when E is along one of the cubic [100] or [010] axes. 
Notice that this is completely different from the spin Hall 
drag current, which is always perpendicular to the elec- 
tric field. However, the spin current can also be made 
perpendicular to E^ by orienting E is along one of the 
principal axes [110] or [110]. 

In the rest of the paper we present some details of our 
theoretical analysis, and provide a numerical estimate 
for the size of the effect. The salient conclusions are as 
follows: (i) The spin transresistivity is proportional to 
{T/EfY which is characteristic of Coulomb drag phe- 
nomena [6l [121 US] (ii) It is inversely proportional to the 



fourth power of the interlayer separation and (iii) It is 
parametrically stronger than the spin Hall drag [6] in the 
clean limit since the transresistivity is proportional to 
Qr, whereas the contribution to the spin-Hall drag is 
independent of r. 

Theory - The part of the "spin-orbit dressed" Coulomb 
interaction which is relevant to the effect described in 
this paper has matrix elements [H [HI [TBj 

M{kiSi,k2S2 -^P1S[,P2S2) = {Vq - yVq((7 + ^5) 
{Xs[Xs^J[^l X (pi^ki)]^ - [&2 X {P2 ^ k2)]z\Xs^Xs2)} 

^fcl+fc2,Pl+P2 5 (I) 

where the subscripts 1 and 2 refer to electrons in layers 1 
and 2, respectively, &i and &2 are the spin operators of 
the first and second electron acting on spinors \Xsi)^Xs[) 
and \Xs2) ^Xs'J ^ respectively. The normalization area is 
set to unity. Here q = p — k is the transferred wave 
vector, 

' ~ ^ {q ^ Qs)^ - q^s^-^^"- 
is the Fourier transform of the screened Coulomb inter- 
action between the layers, L is the distance between the 
layers, z^: is a background dielectric constant and Qs is 
the screening wave vector. The strength of the spin-orbit 
interaction is controlled by the "effective Compton wave- 
length" A for the semiconductor: 

3 i?2(i?^ + A)2' 

where P is the Kane parameter [16 and Eg and A are 
the band-gap and the spin-orbit splitting of the valence 
band. In GaAs 5 A^. Notice that, in writing Eq. ^ 
we have taken into account only the Rashba-like contri- 
bution arising from the component of the Coulomb field 
perpendicular to the layers. The additional term com- 
ing from the component of the Coulomb field parallel to 
the layers was discussed in Refs. \6\ \T7\ and [T8l and is not 
shown in Eq. ([T]). 

As discussed in the introduction and in the caption 
of Fig. [2j a steady current driven in the active layer (2) 
produces, via the "spin-orbit dressed" Coulomb interac- 
tion, a quadrupolar distribution of spin in the passive 
layer. The spin generation rate in the passive layer can 
be calculated by means of the Fermi golden rule with the 
matrix element of interlayer electron-electron interaction 
given by Eq. ([l]) [HIlO]. To first order in A^ this gives 

9k = J, ^k-\-k', p+p'^{£k + ^k' — — £p') X 

^ k'pp' 

Uu-p[{p ^k)xz]{E-{p- k))hfk'{l - fp){l - fp^m 

where Uq = X'^ {q -\- qs)\Vq\'^ ^ E is the electric field acting 
on the carriers in the active layer, m is the effective elec- 
tron mass, ksT is the temperature measured in the en- 
ergy units, Sk = Ti^k'^ /(2m) is the electron dispersion, fk 
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is the Fermi-Dirac distribution function and z is the unit 
vector normal to the layers. We have also assumed that 
equilibrium densities, effective masses and Fermi energies 
are the same in the two layers. In deriving Eq. ([3| we 
took into account only linear-in-£^ nonequilibrium cor- 
rection to the distribution function. It is easy to show 
that, due to the presence of the two factors E • {p — k) 
and + X £, the angular dependence of is an angu- 
lar harmonic of order 2, i.e. we have g^^k sm2(pk and 
9y,k ^ —cos2(pk, where cpk is the angle of k with the x 
axis: this is the quadrupolar pattern of spin generation 
mentioned in the introduction and described in Fig. |2] 

The spin dynamics in the passive layer is governed by 
the spin-orbit splitting of the energy spectrum. Since the 
first harmonic components of the spin splitting (arising 
from linear- in- A;- terms) do not result in a dc spin Hall 
current [TTJ [2TJ [22] , we only take into account the third 
angular harmonics of the Dresselhaus term which is 
inevitable in any zincblende structure. In cubic axes with 
X II [100] and y \\ [010] these have the form 

^x,k = -ftscosScpk, ^y,k = -J^3sin3(/:?fc, (4) 

where = 7c/c^/(2fi), 7c is the bulk Dresselhaus con- 
stant. The electron spin distribution function, s^, is de- 
termined by a kinetic equation, which in the steady state 
takes the form \JJ^ 

Sk 



Sk X 



9k 



(5) 



Here we have assumed that all spin-independent scatter- 
ing processes can be characterized by a single relaxation 
time r. The solution of Eq. ([5| is 

{flk X gk)z 



Sz,k 



1/t2 



(6) 



We note that such a simple form of the solution results 
from keeping only the third angular harmonics in Qk- 
The solution in the general case can be constructed fol- 
lowing Refs. 11 and 23: in such a case the overall spin 
will be smaller due to faster spin relaxation caused by the 
linear-in-/c Dresselhaus term, and the spin distribution 
may have more complicated form due to the anisotropy 
of the spin splitting. 

By definition, the current density of spin z component 

is 



2_^Sz,kVk , 

k 



(7) 



where Vk = fik/m is the electron velocity. Substituting 
the expression for Qk^ obtained from Eq. ([3|, into Eq. (|6| 
for Sz^k^ and then Sz^k in Eq. ([7|), we arrive at our main 
result 

27rfier^ 



kk'pp' 



2Uk-p[^k-{p + k)][E-{p-k)] 
fkfk'{l-fp){l-fp')- 



(8) 



Equation (Isl) can be recast in the standard form: 



(9) 



where Gij is the spin drag conductivity. Our system 
is characterized by the (72v point symmetry group be- 
cause we take into account (i) the Dresselhaus field, 
and (ii) the Rashba-like interaction associated with the 
perpendicular-to-plane component of the Coulomb field. 
In the basis of the principal axes xi \\ [110] and yi \\ [110] 
Eq. (|9| can be written in terms of two independent con- 
stants G and Gi as 

r,^ = {G + G,)Ey,, Jl = {G-Gr)E,,. (10) 



Equation (10) gives the full phenomenological picture of 
the spin drag in the presence of the spin-orbit interaction. 
The drag mechanism described here produces G 7^ and 
Gi = 0. This rather peculiar situation implies that the 
spin current flows parallel to the electric field if the lat- 
ter is applied along the [100] or [010] directions, but, in 
general, it has both a parallel and a perpendicular com- 
ponent. It is only when the electric field is along one of 
the principal axes that we get a pure transverse current, 
but even in this case it is in a sharp contrast with the 
spin Hall drag, namely an electric field along xi produces 
a spin current along yi , but an electric field along yi pro- 
duces a spin current along xi, rather than —xi. These 
unique symmetry-related features will help distinguish- 
ing the predicted new effect. 

We now come to the quantitative evaluation of G. The 
expression ([8| for the spin current density is similar to 
the expression encountered in the calculation of the or- 
dinary drag current [13 . The main difference is that the 
integrand depends not only on the momentum transfer 
q = p — k^ but also on the sum of the initial and final 
momenta, p-\-k. For simplicity we focus here on the case 
of well separated layers, kpL ^ 1, where the momentum 
transfer is typically small {q ~ <C /c^?), and therefore 
one can replace k -\- p ~ 2/ci?, making an error of order 
q/kp. With this approximation, the standard method of 
evaluation can be applied, and, after converting the wave 
vector sums into integrals, and for low enough tempera- 
tures ksT <C Ep^ we obtain 



G 



n^eVt^kFT^ 
' 2T:m?kBT 



J i^^rJ-o. sinh' 



\2kBT ) 



(11) 

where xol'/jw) is the non-interacting density-density re- 
sponse function of the two-dimensional electron gas [24^ . 
Integrating over frequency we get 



G = 



hiqsL) (12) 



with 



hix) = 



SkpL 



f 

Jo 



yHy + x)e-^y 



[{x + yy 



-2y]: 



Ay, 



(13) 
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for A^iT^L > 1. For g^L > 1 the integral Eq. ([Ts]) is 
evaluated to be l2{qsL) ^ 2.96/(/ci7gf L^), as a result 
G (X 1 /L^ as in the ordinary drag effect [HI [13] . 
Results and discussion - It is instructive to estimate the 
drag resistivity which controls the observable spin ac- 
cumulation. This is determined as follows. First, we 
express the external field E via the current density gen- 
erated by this field in the active layer, j = aE^ where 
(J = ne'^r/m is the Drude conductivity of the electrons 
and n is their density. Then, we observe that the spin 
current is associated with some effective electric "spin" 
field E, 



2e 



r = (jE, 



As a result we obtain 



2eJ^ 2eG . 
ha ha"^ 

Hence, the drag resistivity is given by 



2h, 



2 / 2 \ 2 



p. = — (AM^1.3r(M:)'^^) I.fei), (14) 



where vp is the Fermi velocity. Notice that, at vari- 
ance with the (side-jump) spin Hall drag resistivity, the 
present result is proportional to the momentum relax- 
ation time r and therefore it is parametrically dominant 
in the "clean" limit, r ^ oo. 

For the following values of the parameters: electron 
density n = 2 x 10^^ cm~^, r = 40 ps (which corresponds 
to a mobility ja = 10^ cm^/Vs in a GaAs quantum well 
with m = 0.067mo), bulk Dresselhaus splitting constant 
7c = 20 eV A^, = 5 A^, and = 13, we obtain 
(Xkr)'^ - 6 X 10-"^, Qsr - 0.9, and e^/nhvp - 1. Thus, 

taking /2 ~ 1, we have ps ^ ^ i^E^^ Ohm, which for 
(kBT/EF) = 0.1 is 0.04 Ohm, i.e., about two times larger 
than the conservative estimate for the spin Hall drag re- 
sistivity from side jump given in Ref. [6] Such a value 
of spin drag resistivity can be detected by spin Faraday 
or Kerr rotation technique. In InAs based structures the 
spin-orbit coupling parameter is about an an order 
of magnitude higher than in GaAs and Qs is also larger. 
This makes narrow-band semiconductor bi-layers partic- 
ularly suitable for the observation of the spin current 
injection. 

In conclusion, we have described a new coupling mech- 
anism, partly Coulomb and partly spin-orbit, through 
which a spin current can be injected, or a spin accumu- 
lation induced, in an electron layer by a regular electric 
current fiowing in an adjacent layer. The new coupling 
can play a role in the design of circuits in which an elec- 
tric current must be converted into a spin current and 
viceversa. 
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